function [S_min,fval] = UnConstrainedOpt(T,S0)
% 输入目标位置T 和 三个监测站坐标初始猜测值SO
% 输出最佳位置S_min 和 最小GDOP值fval

%% fminunc 无约束多变量函数的最小值
f = @(S)fun_gdop(T,S);
options = optimset('PlotFcns',@optimplotfval);
[S_min,fval]= fminunc(f,S0,options);
% S3_fminsearch = fminsearch(f,S30,options);% fminsearch单纯形法
%% plot 
figure 
plot(T(1),T(2),'*');hold;
plot([S_min(1) S_min(3) S_min(5)],...
[S_min(2) S_min(4) S_min(6)],'o');
% figure
% plot3(T(1),T(2),T(3),'*');hold;
% scatter3([S1(1) S2(1) S3(1)],[S1(2) S2(2) S3(2)],[S1(3) S2(3) 0])
% xlabel('x(km)');ylabel('y(km)');zlabel('z(km)');
% % 平面图
% figure
% plot(T(1),T(2),'*');hold;
% plot([S1(1) S2(1) S3(1)],[S1(2) S2(2) S3(2)],'o');
% xlabel('x(km)');ylabel('y(km)');zlabel('z(km)');